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We describe a detailed study of string-seeded structure formation using high resolution 
numerical simulations in open universes and those with a non-zero cosmological con- 
stant. We provide a semi-analytical model which can reproduce these simulation results 
including the effect from small loops chopped of by the string network. A detailed study 
of cosmic string network properties regarding structure formation is also given, including 
the correlation time, the topological analysis of the source spectrum, the correlation be- 
tween long strings and loops, and the evolution of long-string and loop energy densities. 
For models with V = Qh = 0.1—0.2 and a cold dark matter background, we show that 
the linear density fluctuation power spectrum induced by cosmic strings has both an am- 
plitude at 8h _1 Mpc, o~s, and an overall shape which are consistent within uncertainties 
with those currently inferred from galaxy surveys. The cosmic string scenario with hot 
dark matter requires a strongly scale-dependent bias in order to agree with observations. 

1 Introduction 



One of the outstanding problems in cosmology today is developing a more precise 
understanding of structure formation in the universe, that is, the origin of galaxies 
and other large scale structures. Existing theories for the structure formation of 
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the Universe fall into two categories, based either upon the amplification of quan- 
tum fluctuations in a scalar field during inflation, or upon a symmetry breaking 
phase transition in the early Universe which leads to the formation of topologi- 
cal defects. While techniques for computing density perturbations for the former 
are well established, 1,2,3 only little quantitative work exists for the latter due to 
calculational difficulties in modeling nonlinear effects, especially for cosmic string 
models. 4 

The cosmic string scenario predated inflation as a realistic structure forma- 
tion model, 4 but it has proved computationally much more challenging to make 
robust predictions with which to confront observations. Only until recently, sig- 
nificant progress in understanding cosmic strings as seeds for large-scale struc- 
ture and Cosmic Microwave Background Radiation (CMBR) anisotropics has been 
achieved. 5 ' 6,7 ' 8,9 ' 10 ' 11 ' 12 ' 13 In this paper, we concentrate on local cosmic strings and 
describe a self-consistent method based on a fluid approximation to study string- 
seeded structure formation with either a cold or hot dark matter (CDM or HDM) 
background. The primary quantities of interest for comparing theories to obser- 
vations are the power spectra of fluctuations, both for the mass density and the 
CMBR. In the work presented here, we concentrate only on the former but we will 
briefly discuss the latter. The present work relies on high-resolution numerical sim- 
ulations of a cosmic string network, 14 with a dynamic range extending from well 
before the matter-radiation transition through to deep into the matter era. Sev- 
eral important properties of local cosmic strings are revealed. The resulting power 
spectrum of linear density perturbations, V(k), and the mass fluctuation ampli- 
tude at 8/i _1 Mpc, as, are calculated. A semi-analytic model is also introduced. 
It can reproduce the V(k) of our high-resolution simulations. We further investi- 
gate the dependence of V(k) on the curvature K and the cosmological constant A. 
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In particular, for models with T = flh = 0.1-0.2 and a CDM background, wc 
show that both cr 8 and the overall shape of V{k) are consistent within uncertainties 
with those currently inferred from galaxy surveys. The HDM scenario with cosmic 
strings seems to require a large scale-dependent biasing in order to be consistent 
with observations. 

The framework developed here is also suitable for investigating any other matter 
sources behaving stiffly, whose evolution is largely independent of the background 
matter and any inhomogeneities in the universe. This direct numerical approach 
marks a considerable quantitative advance by incorporating important aspects of 
the relevant physics not included in previous treatments. As such, the cosmic string 
power spectra presented here should be the most reliable to date. 

The structure of the paper is as follows: In section 2, we investigate the linear 
perturbation equations for cosmic defect models in a flat A-universe. We then argue 
that this can be extended to closed and open universes for cosmic defect models. 
Section 3 starts by investigating the topology of generic sources. Then several 
important string properties are described, including the slow relaxation of the long- 
string and loop energy densities from the radiation to the matter era, the topology 
of string network, the string network correlation time, and the correlation between 
long strings and loops, etc. In section 4, we describe the approximation schemes 
invoked in this work, including the compensation of the source into the background, 
the inclusion of HDM, a semi-analytical model which can accurately reproduce the 
simulation results, and a simple extrapolation scheme which can generalize the 
simulation results from a flat A = model to open and flat A-models with any 
desired Hubble parameter. This simple scheme is then numerically verified with 
high accuracy. The unimportance of the late-time non-scaling behavior of strings in 
open or A-models is addressed. In section 5, we present the main results. Empirical 
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formulae describing our main results are presented here. A brief discussion regarding 
the CMBR anisotropies is also given. Finally, a conclusion is given in section 6. We 
have defined the conventions of the background cosmology and given the analytical 
solutions for open, flat, close and A-modcls in Appendix A. Appendix B defines 
the conventions of our power spectrum and variance calculations. Through out the 
paper, we will use h — 0.7, where the Hubble parameter is defined in the usual way 
H = 100/ikms _1 Mpc _1 . A generalization of results from this choice to any other h 
is provided in section 4.4. 

2 Perturbation Equations 

We consider density perturbations about a flat FRW model which are causally 
sourccd by an evolving external-source network with energy-momentum tensor Q , / g(x, 
In the synchronous gauge with a cosmological constant, the linear evolution of the 
radiation and CDM perturbations, S T and S c respectively, are given by (modified 
from Ref. 15 ) 

Sc + -Sc-^(-) 2 (n c S c + 20A) = 47TG6+, (1) 

Ob A \ CI / 

<5r - \v 2 8 r - p c = 0, (2) 

where a(rj) is the scale factor, + = 6oo + 0»», a dot represents a derivative with 
respect to the conformal time 77, and 

" c= l + a + Ba^ + Ca^ ^ = l + a + Ba* + Ca*' (3) 
with B and C defined in (A. 3). Eqns. (1) and (2) are exact when K = 0, i.e. 
Qko = 0. However, they are still a very good approximation even if K 7^ because 
most perturbations on scales of interest are generated when the curvature-associated 
perturbation terms are small. 
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It proves useful to split these linear perturbations into initial (I) and subsequent 
(S) parts, 15 

Mx,rj) =&(x,»j)+$jv(x,fO, N = c,i, (4) 
with initial conditions 

Swim) = Sn(vO, $N(Vi) = Sn(vO, ( 5 ) 

4W) = <^rW) = 0, withr? > jj > 77i . (6) 

In the case of topological defects, the initial perturbations S l (x,i]) depend on the 
defect configuration at the defect formation time 771, because ultimately the forma- 
tion of defects creates under-densities in the initially homogeneous background out 
of which they are carved. The subsequent perturbations 5 (x, rj) are those actively 
generated by the defects themselves for fj > rji. Because defects induce iso-curvature 
perturbations, <5 : (x, 77) must compensate <5 s (x, 77) on comoving scales |x — x'| > 77 
to prevent acausal fluctuation growth on super-horizon scales, as we shall discuss 
later. 

The system of eqns. (1) and (2) can be solved for the subsequent perturbations 
5 s (x, rj) by using a discretized version of the integral equation with Green functions: 

4(x,ry) =4kG J\ J d 3 x'g N (X-,r),r))e + {x',r)), N = c,v, (7) 

where X = |x — x'|. The easiest method of obtaining the Green function solutions 
is to go to Fourier space and solve the resulting homogeneous system of ordinary 
differential equations numerically with appropriate initial conditions. Since the 
Green functions depend only on the modulus of X = |x — x'|, it follows that their 
Fourier amplitudes must only depend on the modulus of k. We will use a tilde ~ to 
denote the Fourier transform of a function. With the change of variable y = 1 + 
(V2 — l)rj/r] e q and using (7), in Fourier space eqns. (1) and (2) with K = A = 
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become 

120r/fi 



(i- y *)g'>-2yg' c + 



y 2 



0, (8) 



4 ~ 4fc 2 ~ 
^-3^ + ^ = 0, (9) 

with A defined by (A. 3), and initial conditions at 77 = 17 (or y = y(fj)) satisfying 

g r = g c = o, g> = p> = ±. (10) 

Here, a prime represents a derivative with respect to y (though only in these equa- 
tions). 

Under certain limits, (8) and (9) can be solved analytically. When kt] » 1, the 
radiation component will oscillate many times per expansion time and will have 
little net effect on the matter. So we can set g r = in this case and thus (8) 
decays to a homogeneous associated Legendre equation, whose solution is the linear 
combination of the two associated Legendre polynomials Pi and Q 2 - When k-q <C 1, 
by (9) and (10) we know £7 r /£> c = 4/3. So (8) decays to a homogeneous associated 
Legendre equation, whose solution is the linear combination of the two associated 
Legendre polynomials P 2 ~ and Q\. One can therefore solve (8) and (9) with (10) 
for g c under these limits to get: 

G (k-n n) — i ffl ~ ^ A ^ B ^ B (v) A ^ > kr ) » 1 > fill 
M i(y*-l)-i[C(y)V(y)-V(y)C(y)}, 1 ' 



where 



B(y) = %?f±, (13) 

c (y) = ^r=- v ( 14 ) 
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We notice that the scale factor 0(77) = y 2 — 1. Because we are only interested in 
the matter perturbations at late times (i.e. today 770) and we know from (8) that 
Q c oc a/a cq = a when r]/r] cq — > 00, what we actually need is the function: 

T c {k;f,) = lim -g c {k;r],fj). (16) 

»7/»)eq— »°0 

Hence for fcr? — > and 00, we obtain from (11) that 

foc(yW)) = |(y 2 -lM(y), fcr ? »l, (17) 
T (y(fi)) = ^C(y), fcr/«l. (18) 

Fig. 1 shows how (17) and (18) evolve with time. Fig. 2 shows the numerical so- 




-i -P _1 n 1 p 
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Fig. 1. The evolution of the functions T c (r/) = T c (fc; j)) under the two limits: kr\ <C 1 (dashed line), 
and fc?j S> 1 (solid line). 

lutions of Q c (k;r]o,fj). It confirms the asymptotic behaviors in these two regimes: 
on super-horizon scales (k -C 277/77), Gc(k;rjo,7j) scales as rj -1 , which is indicated 
by (18); on sub-horizon scales (fc ^> 2n/t]), the growth of matter perturbations has 
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a maximum at r\ ~ r] cq , as seen in (17) (see also Fig. 1). We note that on sub- 
horizon scales, when going from 77 = 0.057i cq to f\ = 5r7 cq , G c {k;r]o,fi) changes only 
within a factor of 2. Adding the fact that cosmic defects seed matter perturbations 
only on sub-horizon modes because on super-horizon scales they are compensated 
by the background and therefore can not create density perturbations, we see from 
Fig. 1 and 2 that defects induce power into matter perturbations mainly during 
the radiation-matter transition regime. This is generically different from inflation- 
ary models, in which matter perturbations are seeded during inflation in the deep 
radiation regime when all the modes are well outside the horizon. 

Now we can obtain an approximated solution for the Green function by combin- 
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ing these two modes: 



Gc(k;r]o,fl) 



a T c (k;fj) 



a {Tco + (T - Too) [1 + (afcW 4 / 3 ] '} 



with 



( a = 0.75, for 0.03 < rj/r] eq < 125 and 3/iMpc" 1 < k, 
\ a = 0.5, for 0.2 < r)/ry cq < 125. 



(19) 



The combination factor in (19) is a numerically verified guess and its accuracy is 
plotted in Fig. 3. We note that the late-time positive departure of the accuracy (e.g. 
fj = 100, 125?7oq, see Fig. 3) is caused by neglecting the negative terms associated 
with B(y) and V{y) in (11) when taking the limit (16) to get (17) and (18). 

Finally, we obtain the subsequent perturbations in Fourier space today by in- 
tegrating the Green function with the external source throughout the dynamical 
range within which the external source is present: 



Here the Green function G c {k;r)o,fj) can be either (19), or numerically obtained 
from eqns. (8) and (9) with initial conditions (10). 

We also tested the performance of the fluid approximation (1) and (2), by com- 
puting the "transfer function" for primordial density perturbations in the absence 
of a source (Tfc), defined as the ratio of the growing mode coefficient deep in the 
matter era to that deep in the radiation era. In Fig. 4, we compare our result with 
other previous work for the case of a flat universe with no cosmological constant. 
We see that although we have used the fluid approximation to model only the mat- 
ter and radiation components, our result is only slightly larger with a maximum 
deviation being 5% from the fit given by Bond and Efstathiou (BE), 17 and 8% from 
the result from CMBFAST, 18 on the scale k ~ 0.5/i 2 Mpc _1 . We have also verified 
that the performance of the fluid approximation is very weakly dependent on the 




(20) 
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choice of cosmological parameters. 
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Fig. 4. The transfer function of standard CDM model. The solid line is our result using the fluid 
approximation (1) and (2); the dashed line is from CMBFAST, 18 with Q c0 = 0.97, Q b = 0.03, 
and h = 0.75; the crosses show the fit by Bond and Efstathiou (BE), 17 with the same parameter 
choice. 



3 The Perturbation Source — Cosmic Strings 



3.1 Topological analysis 



The types of external sources we will come across can be roughly classified into 
four categories, according to their topology: zero-, one-, two-, or three-dimensional 
objects. For the first three types, Ref. 19 gave the forms of their power spectra by 
Fourier transforming a set of randomly distributed points, straight filaments, and 
spherical shells: 



V 2 



-Pi(k) 



V 2 tan-iQc) 
8tt 3 kc* 



(21) 
(22) 
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V ^ k) = &<xrm®' (23) 

where V is the volume of the box to be Fourier transformed, and we have taken the 
distributions p(c) = c 2 exp(— c/c*)/2c^ for filament lengths &ndp(R) = R 2 exp(— R/R if )/2Rl 
for shell radii. This leads to (c) = 2\/3c* and (R) = 2\/3-R*. Here, we further 
consider three-dimensional objects. The power spectrum of a set of randomly dis- 
tributed uniform spheres, with a radius distribution ofp(r) = r 6 exp(— r/r»)/720r£ 
(which gives an rms radius (r) = 2\/l4r*) is: 

V 2 (l + 0.8fc 2 r 2 ) 
W = ^3 (1 + 4^)3 • ( 24 ) 

We can see that for sufficiently small k, they all give a white-noise power spectrum 
V(k) — V 2 /8ir 3 oc fc°. This is because on scales sufficiently above their characteristic 
scales 2v / 3c», 2\/3-R* and 2\/l4r*, all objects are essentially point-like. On the 
other hand, at larger k, they give the characteristic behaviors of zero-, one-, two- 
and three-dimensional objects: V(k) oc k°, k^ 1 , k~ 2 and fc -4 respectively. In Fig. 5, 
we show the power spectra of one realization of a set of randomly distributed zero-, 
one-, two-, and three-dimensional objects. They agree very well with the analytical 
behaviors of (21), (22), (23) and (24). The reason for having a steeper small-scale 
slope in the power spectrum of higher-dimensional objects with the same overall 
mean density is that the mass is more diluted in this case. 

3.2 Cosmic strings 

The external source under consideration is a network of cosmic strings. For local 
gauge strings, the microphysical width of the string will be many orders of magni- 
tude smaller than its typical curvature radius, so that we can take the zero thick- 
ness limit. Therefore, the cosmic strings which dynamically source the subsequent 
perturbations 5 s (x, rf) in (4) have spacetime trajectories which we can represent 
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Q_ 




Fig. 5. The power spectra of randomly distributed zero, one, two and three dimensional objects 
(dashed lines, downwards). Their characteristic sizes have been normalized to 2n, and we have 
also used (27) to deconvolve the lattice effect. These spectra agree very well on small scales with 
those solid lines, which have exact slopes of 0, —1, —2 and —4 respectively. On large scales, they 
all give white noise V(k) <x k°. 
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as x£ = (rj, x s (cr, 7])), where a is a space- like parameter labeling points along the 
string. The stress energy tensor of the string source is then given by 4 

9 Q/3 (x, »?)=/*/ da(ex«xP - e^x'^x'f )<5 3 (x - x s (a, rfi), (25) 

where /i is the string linear energy density a prime represents a derivative with 
respect to a, e = [x' s 2 /(l — x 2 )] 1 / 2 , and we have also taken that x s • x s ' = 0. It is 
then straightforward to compute 6+ in (1) as 

e+(x, n) = e 00 + e« = 2^ j dcrex 2 <5 3 (x - x s (a, ??)). (26) 

This source term can be calculated directly from the string network which was 
evolved using the Allen- Shellard (AS) string simulation, 1 and it is in good agree- 
ment with the other high resolution simulation of Bennett-Bouchet. 20 ' 21 Large-scale 
parallelized simulations were performed on the COSMOS supercomputer, a Silicon 
Graphics Origin2000 with 20 Gbytes main memory. Dynamic ranges exceeding one 
thousand expansion times were feasible because of the simulation size and by using a 
modified 'point-joining' algorithm which maintained fixed comoving resolution. Be- 
cause sampling points are used to trace cosmic strings, we invoke the "cloud in cell" 
(CIC) method to assign these sampling points of strings onto the comoving grids 
and then Fast Fourier Transform them. We then multiply this transformed quan- 
tity on grids by the following sharpening function to compensate for the smoothing 
effect induced by the CIC method: 22 



(fc/2) 2 

sin 2 (fci/2) + sin 2 (fc 2 /2) + sin 2 (fc 3 /2) 



^) = .,2 t , I • 2 , ^ I ■ 2n- ^ > ( 27 ) 



where k = (ki, &2, Aft) and k = |k|. This scheme has been verified against sets of 
zero, one, two and three dimensional objects as we have seen in Fig. 5. 

We will divide the string network into two parts: the long strings and small 
loops. We define the latter by taking those close loops whose lengths are smaller 
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than one tenth of the horizon size, and we find that the resulting energy evolution 
in either the long strings or the loops is not very sensitive to this threshold because 
the loop sizes are typically at least two orders below the horizon size and there are 
not many large loops around the threshold. We will first investigate the properties 
of long-string network in this subsection, and leave the investigation of loops to the 
next subsection. 

Fig. 6 shows the initial and scaling long-string configurations, and their power 
spectra. Because the initial network is a random walk with a correlation length £ 
along the x, y or z directions, its properties are very similar to that of a sampling 
function (the periodic Dirac Delta) 

oo 

II(aO= 8(x-n0, (28) 

n— — oo 

where n is an integer. We know its Fourier transform is also a sampling function. 
Hence, the power spectrum of this string network has periodic peaks, whose interval 
length relates to the correlation length as Ak — 2tt/£_. Due to the 3-dimensional 
structure, there are also other minor peaks corresponding to 2\/2n/£, 2\/37r/£, etc. 
With £ = 77/8, we also note that at small scales (k > 1677/77), the power spectra 
have an overall slope of —1 because strings are line-like objects; at large scales 
(k < 1671-/77), the power spectra turn flat because of the point- like property (sec 
subsection 3.1 for the topological analysis). 

On the other hand, the scaling string network in Fig. 6 shows different properties. 
Its power spectrum is much smoother than the previous case because the periodicity 
of strings has no preferred directions (i.e. the x, y and z directions in the previous 
case) and therefore the periodic peaks degenerate to only one broad peak here, 
corresponding to the mean spacing of strings. This smoothness feature of the power 
spectrum can be employed to ensure that when the simulation of structure formation 
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log 10 kr| 

Fig. 6. The top-left is an initial network of cosmic strings, the top-right a scaling long-string 
network, and the bottom their power spectra. The first peak in the initial spectrum appears at 
fc^ = Wn/n, which accurately corresponds to the initial comoving correlation length £com = n/8 
set by hand. The turnover scale in the scaling regime is clearly fcj = 20/r/, and the slope of the 
spectrum turns from —1 in the initial network to —1.13 in the scaling regime. The normalization 
here is arbitrary. 
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starts, the string network is already scaling so that the initial configuration of the 
network will have no effect on the subsequent result. We also note that instead 
of turning from k° to fc _1 very smoothly at the correlation scale as we have seen 
from randomly distributed filaments (see Fig. 5), the power spectrum has a broad 
peak on that scale. This illustrates the fact that in the scaling regime, a real string 
network is not like random walk, but has a configuration in which all strings remains 
a roughly constant distance (which is about the correlation length) from each other, 
so that the power on this scale is amplified to form a broad peak. We verify that this 
broad peak appears at k^ w 2n/^ rts 20/ry persistently from deep in the radiation 
era through to deep into the matter era, with an amplitude of about three times 
larger than the large-scale white noise. Here £ w 77/3 is the comoving mean distance 
among long strings. Therefore, any previous work based on modeling cosmic strings 
with random filaments would have underestimated the resulting power spectrum by 
a maximum factor of three (e.g. Ref. 10,23 ), because such a broad peak in a real 
string network appears persistently with an amplitude of about three times larger 
than the large-scale white noise. 

We further notice in Fig. 6 that the slope of the spectrum in the scaling regime 
is not exactly —1 on small scales, but —1.13. This is because strings are very wiggly 
on small scales and therefore have an overall dimension of slightly greater than one 
so as to steepen the small-scale slope to —1.13. 

Another important property of the string network is the correlation time. Be- 
cause the strings' configurations are uncorrelated when separated by a sufficiently 
long time, we can define a correlation time f] c {k, 77) as the time period above which 
the strings' configurations are uncorrelated, i.e. the unequal time correlator (UETC) 



(6+(k, 77)6+ (k, 7/)) w for |r? - rf\ > Vc ■ 



(29) 
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We see that on smaller scales inside the horizon the correlation time is shorter, and 
that on the horizon scale the correlation time r] c is about the order of the conformal 
time rj, although the UETC is dominated within a dynamic range rf = (1 ± 0.2)77, 
taking the half-maximum threshold. In a quantitative analysis, an accurate cal- 
culation of UETC's requires a dynamic range of at least 4 in conformal time. 
This is equivalent to a dynamic range of 16 in radiation-era physical time, and 
64 in matter-era physical time. So far in the literature, only our simulations have 
achieved this criterion for local cosmic strings, while keeping a highest resolution 
equivalent to 1000 3 points in a comoving box. Each of our simulations typically has 
a dynamic range of more than 10 in conformal time. We also find in the UETC 
(0+(k, 7i)8 + (k, rf)) that a hot spot locates at around rf = r\ and k ss 20/ri, 24 which 
is exactly the correlation wavenumber k^ we have previously seen in the string source 
power spectrum. If we take a slice through rf = rj in the UETC (0+(k, ?y)0+(k, 77')), 
then the profile of this slice gives the scaling spectrum in Fig. 6, with a broad peak 
corresponding to the hot spot in the UETC. 

Fig. 7 shows the integration of the UETC over time. This quantity reflects the 
topology of the path swept out by strings. Since a two-dimensional object will have 
a small-scale slope of —2, we expect such a quantity to behave in the same way. 
However, as we can see from the figure, the small-scale slope is not exactly —2, but 
—2.25. This is due to the string wiggliness which is reflected in the wiggliness of 
the wakes they generate. We investigated the level of such a wiggliness in wakes, 
and verified that the wiggles of cosmic string wakes are typically about one tenth 
of the correlation length. 

Another important aspect regarding the cosmic string evolution is the radiation- 
matter transition process. The key observation was the very slow relaxation to the 
matter era scaling density. Fig. 8 shows the evolution of long-string energy density 
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Fig. 7. The integrated UETC, f(0+(k, J?)©+(k, 7]'))dvi' (solid line), and the power spectrum of 
wakes with a thickness of 1/10 of their correlation length (dashed line). The normalization is 
arbitrary. The dotted and dot-dashed lines have exact slopes of —2 and —2.25 respectively. 
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Poo. A good analytic fit to the evolution of ((rj) = poet 2 / p is found to be 




C(??) = 14.7 : pr 

When compared with simulations of higher resolutions, we find that our values here 
are about 10% higher, while the overall shape remains unchanged. This is due to a 
lower kink density in our case, shortening the mean length of strings. We can see 
from Fig. 8 that the relaxation process extended well beyond recombination (for 
h = 0.7, Iog 10 (r7 roc /f7 oq ) w 0.76 and log 10 (?7o/r7cq) ~ 2.4). This has important impli- 
cations for large-scale structure and CMBR anisotropics as we shall discuss. Fig. 8 
also shows that the long string energy density at the radiation-matter transition is 
closer to the radiation era asymptotic value than previously thought. 



(30) 
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3.3 Cosmic string loops 

The role of small loops produced by the string network has evolved from a poten- 
tial one-to-one correspondence between loops and cosmological objects, 25 through 
to a completely subsidiary role relative to the wakes swept out by long strings. 26 
This dethronement of loops was a result of numerical studies which showed that 
the average loop size I — at was much smaller than the horizon, I « d^; 14,27 
they may even be as small as the lengthscale set by gravitational back-reaction 
a ~ 10~ 4 , a value appropriate for GUT-scale strings. 4 Add the high ballistic loop 
velocities observed v sa and it was not surprising that these tiny loops have 

been assumed to be more or less uniformly distributed and hence a negligible source 
relative to the long-string network. 19 Nevertheless, small loops always make up a 
significant fraction of the total string energy density at any one time and, as we will 
demonstrate later, loop-induced inhomogencitics arc considerable if their lifetime 
is not much smaller than the Hubble time. By properly incorporating these loop 
perturbations, we will show that their contribution relative to the long-string wakes 
is almost comparable and also highly correlated with these wakes. 

To investigate loops, we keep all the loops generated by the long-string network, 
and model those which are smaller than a fixed fraction of the horizon size as rela- 
tivistic point masses. The effects of the evaporation of these loops into gravitational 
waves and the damping of loop motion due to expansion are also included. Note 
that here we assume the standard scenario with an evolving 'infinite' string network, 
but there are alternative models in which the initial string configuration consists 
entirely of loops (e.g., Ref. 28 ). 

To begin with, let us consider some analytical properties of loops. First, the 
Nambu equations of motion for cosmic strings in an expanding universe can be 
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averaged to yield: 

+ 2H(l + ( V 2 )) Poo = -X L , (31) 



dt 

where (v 2 ) is the mean square velocity of strings and Xl is the transfer rate of 
energy density from long strings into loops. In the scaling regime the long-string 
energy density should scale with the background energy density evolving as 

dpoo = _ 2 Poo (32) 
dt t y ' 

Substituting this into (31) to eliminate dpoo/dt gives 

tXL _ ( (1 — (v 2 }) ~ 0.6 in radiation era. 



§(l-2«))~0.2 in matter era, ^ 

where (v 2 ) > (v^) <~ 0.6. 14,27 Both (32) and (33) provide a check for the scaling 
behavior of long strings and loops in the cosmic string network simulations. Figure 
9 shows the evolution of Xl- We can see that the expected amount of energy 
was converted into loops in our simulations so that X L has the correct asymptotic 
behavior given by (33). However, the typical loop-size (and consequently their 
lifetime) does not approach scaling so rapidly and is therefore larger than physically 
expected for most of the duration in the simulations. 14,27 To overcome this problem 
we need to rescale the loop lifetime in an appropriate way. 

The loops produced by a cosmic string network will decay into gravitational 
radiation, with a roughly constant decay rate TG/j, 2 , where fi is the string linear 
energy density. Typically T = 50 — 100 with an average (r) <~ 65. 29,30 In very high- 
resolution simulations, we observe a loop-production scale on which the loops are 
most effectively generated by the network, although this is only approached towards 
the end of simulations used here. For both analytical and numerical simplicity, we 
make the reasonable assumption that the loop production as 'monochromatic', that 
is, all loops formed at the same time will have the same mass in the limiting scale- 
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Fig. 9. Evolution of X^(t). The dotted lines are the asymptotic values in the radiation and the 
matter eras in (33). 
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invariant distribution. We can write the initial rest mass of a loop formed at time 

M* L = aiiU = fYGn 2 U ■ (34) 

Here, the parameter / = a/TGn is expected to be of order unity if the size of the 
loops formed at the time t is determined by gravitational radiation back-reaction, 
which smoothes strings on scales smaller than TG^it. Thus the uncertainty in the 
average mass and therefore the lifetime of loops formed at a given time is quantified 
by the choice of the parameter /. We also find that the initial rms velocity of loops 
observed from the simulations is (v 2 ) 1 / 2 > 0.7c throughout all the regimes. 

In the simulations we only impose the monochromatic assumption through / on 
the loop decay lifetimes. Thus we have the rest mass of a loop formed at time t* 
evolving as 

M L (U,t) = M* L W{U,t) , (35) 

where 

I otherwise. 
Here r(£*) w /t» is the lifetime of loops produced at time t* (/ = 2, 3 implies the 
decay occurs in one horizon time in the radiation and matter eras respectively). 
The evolution of the loop energy density is then given by: 



PL(t) 



Jo L a{t) 



W(t',t)dt' 



*(*) 

f/t 2 for / « 1, 

{ Vf/t 2 for /> 1 (radiation era), (37) 
(ln/)/t 2 for / > 1 (matter era), 

where we have used the scaling behavior (32) and (33). Consequently, the scaling 

of the power spectrum induced by loops in / should interpolate between / 2 and / 

(radiation era) or (In/) 2 (matter era). We notice in (37) that we have ignored the 

effect of loop velocity redshifting due to the expansion of the Universe, which causes 
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a change in the effective mass. Because loops are formed with relativistic velocities, 
we expect this damping mechanism to have the strongest effect for / 3> 1, but to 
be negligible for /< 1. 

If a loop formed at time t* has an initial physical velocity v», its trajectory in 
physical space accounting for the expansion of the Universe is then given by: 

x(i) = x(t.) +a(t) [ , A =-^7T , (38) 

for t > i*, where A = 7*v*a*, A = |A| and 7* = (1 — Iv*! 2 ) 1 / 2 . Here we have 
neglected the acceleration of loops due to the momentum carried away by the grav- 
itational radiation, the so-called 'rocket effect'. A numerical calculation for several 
asymmetric loops shows that the rate of momentum radiation from an oscillating 
loop is 

\P\=T P G^, (39) 

where Tp ~ 10. 31 Combining with (35), one can show that this rocket effect will 
become important only when: 

l^ + TTT77WWr < 40 > 

which affects only the final stages of the loop lifetime as long as Tp/(Tjv) < 1, 
or equivalently v > 0.15c. For a typical ~ c/\/2, one requires a loop lifetime 
> 43t* in the radiation era and > 17i* in the matter era to rcdshift down to this 
critical velocity according to (38) . Since the values of / we explore here are of order 
unity, it is a reasonable approximation to neglect the transfer of momentum due to 
gravitational radiation. 

With the treatment of (34) and the effects of (35) and (38) Fig. 10 shows a 
loop spectrum with / = 1 in the scaling regime. As we can see, on small scales 
(krj < 100) the shape of loop spectrum is identical to that of long strings (see 
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0.5 1 1.5 2 2.5 

log 1Q kr| 

Fig. 10. On the top is a snap shot of loop distribution in the scaling regime. The size of each point 
reflects its mass. The box size is the same as the horizon. On the bottom is their power spectrum 
with arbitrary normalization. The dot-dashed lines have exact slopes of —1.13 and —0.21. The 
dotted line indicates the correlation wavenumbcr = 20/n seen in the long-string network. On 
large scales (krj < 100), we see loops mimic long strings so that their spectrum shape is identical 
to that of long strings (see Fig. 6). The normalization is arbitrary. 
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Fig. 6). This is because loops trace the paths of long strings even after they are 
formed. We verify this by calculating the correlation coefficient between long strings 
and loops in the scaling regime. This coefficient is always greater than 0.5 on larger 
scales. This phenomenon has also been verified by observing the movies of string 
evolution made from the simulations. On the other hand, on small scales the loop 
spectrum is not exactly a white noise. It has a slope of —0.21. This is because 
even on small scales loops are still not exactly point-like, but with some clumpy 
structures, so that their spectrum slope is steepened as discussed in subsection 3.1. 
In section 5.4, we will see further that this correlation between the long-string and 
loop distributions can have a significant effect on the resulting power spectrum of 
matter perturbations. 

4 Approximation Schemes 
4.1 Compensation scale 

It is a very substantial numerical challenge to evolve the initial and subsequent 
perturbations induced by cosmic strings in (4) such that they accurately cancel on 
super-horizon scales by the present day rj . Long-string network scaling entails the 
copious production of enormous numbers of tiny loops, 14,20,21 whose evolving dis- 
tribution and decay must be carefully followed as described in subsection 3.3. For 
the large dynamic range required for the present study, however, we have by neces- 
sity adopted the compensation factor approximation suggested in a semi-analytic 
context in Rcf. 32 . To implement this, we accurately evolved the string network nu- 
merically, and then multiplied the Fourier transform of the resulting stress energy 
0+(k, rf) by a cut-off function F(k,rj) given by 




(41) 
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where k c (r]) is the compensation wave number. This results in the correct k 4 fall-off 
in the power spectrum at large wavelengths above the compensation scale k~ x ~ 77. 
Thus we obtain 

j c (k, m) = sl(K no) + s s c (k, no) « f g c (t, %, »/)e + (k, v')F(k, v ')d v ' . (42) 

J Hi 

Several proposals for k c (rj) have been discussed in the literature. In Ref. 32 , 
k c = 2nrj~ 1 was suggested as physically plausible, but not seriously justified. More 
recently, however, in Ref. 33 the efficacy of the approximation (42) has been demon- 
strated by studying multi-fluid compensation back-reaction effects in greater detail. 
It is claimed that the compensation scale arises naturally and uniquely from an alge- 
braic identity in the perturbation analysis. For the present study we have adopted 
the analytic fit for k c (r]) presented in Ref. 33 : 

_ V6(3AV + I2^+16) 

where A is defined by (A. 3). This k c (rj) smoothly interpolates from k c = 
in the radiation era to k c = ^/lSr/^ 1 in the matter era. Of course, there is some 
dependence of the perturbation amplitude on this choice. Fig. 11 shows the effect 
of different compensation scales on the resulting power spectrum. Here we have 
used the semi-analytical model which will be described in section 4.3 and verified in 
section 5.1. Evidently, a larger compensation scale (at the same conformal time) will 
give significantly more power at larger scales. For example, with the smaller matter 
era k c of Ref. 32 , the amplitude is approximately 40% smaller around length-scales 
of 60/i _1 Mpc and larger. 

Clearly, the quantitative implementation of compensation effects is one of the 
key uncertainties in all previously published work on gauged cosmic strings. The 
new analytic compensation factor approximation which we use here should improve 
this situation, but the problem ultimately requires a full-scale numerical treatment 
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Fig. 11. A comparison of spectra using different compensation scales. The solid line used the 
definition in Ref. 33 , where \c varies from \/6 in the radiation epoch to \/l8 in the matter epoch; 
the dashed line used a compensation scale suggested in Rcf. 32 . The dynamical range is r) = 
(0.01,50)r, oq . 
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in which all strings and background fluid components are accurately evolved through 
to the present day. 34 



4.2 Hot dark matter 

In order to study the formation of structures with cosmic strings in HDM models 
we use a straightforward modification to the perturbation source similar to that 
employed in Ref. 35 . This method is a reasonably accurate alternative to much 
more elaborate calculations using the collisionlcss Boltzmann equation with de- 
fect sources. We simply multiply the Fourier transform of the string source term 
Q + (k,rf) by a damping factor G{k,rj) given by 



G(k, v ) 



-l 4.43 

1 



1 + (0.435fcD(?7)) 203 



(44) 



where D{r}) is the comoving distance traveled by a neutrino with momentum T v /m v 
from time r\ to 00. The factor G(k, 77) is a fit to numerical calculations of the transfer 
function of a Fermi-Dirac distribution of non-relativistic neutrinos, and accounts 
for the damping of small-scale perturbations due to neutrino free-streaming. 35 We 
calculated D(r)) numerically and found an excellent fit (within 2% error, see Fig. 12) 
for T„ = 1-6914 x 1CT 13 GeV and m v = 91.5 n h0 h 2 eV: 

£(r?) = ^log(^ + l). (45) 



4.3 A semi-analytical model 

The other key difficulty facing defect simulations is their limited dynamic range. At 
any one time, an evolving string network sources significant power over a length- 
scale range which exceeds an order of magnitude. Hence, even for simulations 
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Fig. 12. Neutrino free-streaming length D(rj) with T„o = 1.6914 X 10~ 13 GcV and rn v = 
91.5 Oho'i 2 eV. The analytic fit (45) agrees with the numerical solution within a maximum of 
2% error. 

with a dynamic range of two orders of magnitude in conformal time, the power 
spectrum will only be reliable over one order, even before taking box-size limita- 
tions into account. Fortunately, however, we can invoke a semi-analytic model to 
compensate for this missing power, 19 ' 32 which proves to be fairly accurate in the 
scaling regimes away from the radiation-matter transition. The procedure is essen- 
tially to square the expression (42) and then average over directions to obtain the 
power spectrum V(k). This becomes a Green function integral over the UETC's 

(e + (M)e + (k',7/)>: 32 

(6+(k, r?)e+(k', rf))F(k, V )F(k', ^drjdrf. (46) 

This can be simplified by noting that although there are two time integrals, the 
only significant contributions come from times when 77 and 7/ are reasonably close. 
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That is because the strings' configurations are uncorrelated when separated by a 
sufficiently long time, which is the correlation time r] c (k, rj) as defined in (29). Thus 
(46) can be simplified as 

fi° „ 
V(k) = l6ir 2 G 2 LL 2 \G c {k;rio,v)\ 2 Hk,v)F 2 (k,v)dv, (47) 

where 

G c (k;Vo,v)F(k 1 rj) 
» /*<e + (k,r,)e+(k^WW- (49) 

J -77c 

The second step in cqn. (49), is a reasonable approximation for the following reasons. 
First, the correlation time rj c scales with rj and is at most rj, while the UETC is 
mainly contributed from a dynamic range of (1 ± 0.2)77 taking the half-maximum 
threshold as mentioned in subsection 3.2. Second, around r^q, Gc{k';rfo,Tf') changes 
only within a factor of 2 when going from 0.05rj cq to 5rj cq on sub-horizon scales. 
The structure function T{k, rj) can then be obtained as discussed below: 

(i) The form of !F(k,rf): Cosmic strings are line- like objects which move rela- 
tivistically and so the trajectories they sweep out will be two-dimensional. 
According to the study in section 3.1, we know sheet-like objects have power 
spectra T(k) oc k" 2 at small scales and T(k) oc k° at large scales. For cosmic 
string wakes, we get V(k) oc k~ 2 25 on small scales because of their wig- 
gliness (see subsection 3.2). According to (48), because both F(k',rj') and 
Gc(k';r]o,ri') are scale-independent on small scales (see Fig. 2), the slope of 
!F{k,rj) on these scales should be identical to that of (49). We have seen this 
slope to be —2.25 in Fig. 7. Although in the radiation-dominated regime, 
strings are more wiggly and therefore have a steeper small-scale slope than 
—2.25, we found it sufficient to use this value throughout for the current pur- 
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poses. In addition, because cosmic strings scale with the horizon size, the 
turnover scale from V(k) oc fc~ 2 - 25 to V(k) oc k° scales with the conformal 
time. Adding all these factors together, we can schematically write down the 
structure function in the form: 

HKV) = £{V) 22S/ = £(V) 22 ,/ , (50) 

where £{rj) is the overall normalization against time, / is a constant, fc to = 
(Ii])^ 1 the turnover wave number, and the power n controls the sharpness 
of the turnover behavior. We notice that the broad peak seen in Fig. 7 is 
not modeled in (50). However, we will see in subsection 5.1 that such a 
simplification has a negligible effect on the final power spectrum once the 
phenomenological structure function !F{k, 77) is accurately calibrated by high- 
resolution simulations over a large dynamic range, 
(ii) The time dependence: A convenient way to investigate the time dependence 
£(-q) is to look at the horizon scale k w r/ -1 , at which !F(k, 77)1^=1 ~ £{v)- 
Scaling implies that: 



= a 2 ap B «aVp£g) 1/2 «'T 1/2 - (51) 

hor 



dr] n 

hor 



(52) 



Thus by (49) the structure function is 

^(M)lhar « J^I-l^+MQ+iKV + V")) 

« !\ V- 1/2 (V + v")- 1/2 dv" = 2 -f + 2 (f ) • 
Since cosmic strings scale with the horizon size in the scaling regime, the ratio 
j] c /t] is a constant and therefore: 

£{rj) w T{k, ?7)|hor ~ constant , (53) 

deep in the radiation or matter era. This result is consistent with that presented 
in Ref. 32 , but here we give a clear physical reason for the behavior of the structure 
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function and obtain an explicit expression for the relationship between !F(k, if) and 
the correlation time. The longer the correlation time, the larger the normalization 
of f(k,r]). Therefore an insufficient dynamic range when calculating the UETC's 
or p(k,rf) will result in a severe under-estimate of the final power spectrum. 

Instead of taking the approximation (49), we use simulations to calibrate F{k, rj). 
This is a more accurate procedure specially in the matter era because G c (k'; 7/0, rf) = 
Gc(vo, v') decays linearly with rj'. By fitting the shape and amplitude of the power 
spectrum calculated from (47) with simulations of limited dynamic range deep in 
the matter and radiation eras, we estimated T(k, 77) in these two regimes. We were 
then able to use an interpolation based on the actual behavior of the string density 
during the transition era (see Fig. 8) to provide the normalization of T(k, 77) which 
interpolates smoothly from the radiation to the matter era. This turns out to fit 
well the shape of the actual evolution of £{rf) in (50). We found / = 1/20 and n = 8 
in (50), and the resulting structure function is: 

2 . (54) 

In (54), we note first that the turnover scale is the same for both the radiation 
and the matter eras, i.e. k to = 20 /i]. This corresponds to exactly the correlation 
length of long strings which we have seen in subsection 3.2. Second, £ (77) in (54) 
interpolates smoothly from 0.028 deep in the radiation era to 0.011 deep in the 
matter era. Now with Q c (k; 770, 77) numerically calculated from eqns. (1) and (2), 
we can use the eqns. (47) and (54) to extrapolate the lost power due to the limited 
dynamical range in our numerical simulation. We notice that we have not considered 
the effect from loops here which will be discussed in subsection 5.4. 



F(k,ri) 



fl + (W20V 



lO. 28125 



with £(rj) = 0.028 



0.37322 



1 + (4??eq/r?) 



1.1 
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4.4 The dependence of S(k) on fi,„o, A, and h 

Since an open universe is strongly favoured by many observations, and the existence 
of a cosmological constant is studied in many inflationary models, 36,37 it is natural 
to explore the spectrum in these two regimes. To do this, we first introduce a simple 
rescaling scheme to extrapolate the simulated f2 m0 = 1 and A = spectrum to open 
and A-models, and then verify its accuracy. The rescaling scheme is (adapted from 
Rcf. 9 ): 

s(k, h, o m0 , r» A0 ) = s(k, i, i, o) x n 2 ma h 4 x / 2 (Q m0 , n A0 ) x ff 2 (n m0 , n A0 ), (55) 

where k is in units of O m0 /j 2 Mpc -1 , /(f2 m o,f2 A o) and g(n m0j n A0 ) are given by 

f(0 \ - / fi m0 3 for O m0 < 1, f^AO = . , 

5(^ m o,0 AO ) = ^ =-. (57) 

2 [o^ 7 - n A0 + (l + n m0 /2)(i + n A0 /70) 

In (55), the leading factor fl^h 4 results from the fact that the ratio of scale factors 
a /a cq is proportional to Sl m0 /i 2 and that the Green function (19) is proportional 
to this ratio. The middle factor /(17 m0 ,^ A o) reflects the dependence of the COBE 
normalization of G\i on the cosmological parameters n m0 and f2 A o- It changes 
upwards as we decrease the matter density in an open or A-universe. The last 
factor g(QmO, ^ao) takes into account the fact that in an open or A-universe the 
linear growth of density perturbations is suppressed relative to an f2 m0 = 1 and 
^ao = universe 38 (also verified in Ref. 39 for primordial perturbations). There is 
also a rescaling of k implicit in (55) (the unit of k here is in O m0 /i 2 Mpc -1 rather 
than /iMpc -1 ). This rescaling is due to the fact that the horizon size at radiation- 
matter equality is proportional to (fimofe 2 )' 1 . Hence the physical grid spacing in 
our simulation should be rescaled by this factor, since our length scale is in the unit 

of 7] cq . 
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To verify the accuracy of this scheme, we use (47) together with (54) and the 
Green functions obtained from solving (1) and (2) numerically, to get the spectrum 
S'num(fc) for various choices of fi m0 and Q.ao- We then compare this with S ap (k) 
which is extrapolated from a K = A = model using (55). The initial time of the 
integral is set to r}\ — 0.1, at which both the curvature and A terms are negligible, 
and the final time is today. Fig. 13 shows that on the scales of interest, the rescaling 
scheme (55) is accurate within few percent error for reasonable choices of f2 m o and 
£1ao- A similar independent work regarding the accuracy of (55) in generic defect 
models is also done in Ref. 40 . 

Fig. 14 shows the relative contribution of different epochs to the total power 
spectrum. We can see that most power on the scales of interest is seeded around 
77eq, especially for models with fi m0 < 1. As illustrated in section 2, this is because 
at early times in the radiation-dominated regime, the growth of small-scale pertur- 
bations are suppressed by pressure, while at late times in the matter-dominated 
regime, the large-scale perturbations have less time to grow. This implies that the 
late time non-scaling behavior of cosmic strings has a negligible effect on matter 
density perturbations on scales of interest. We verified that for wavenumbers larger 
than k ~ 0.01/iMpc -1 the late time non-scaling strings contribute less than 3% to 
the total power spectrum for any reasonable choices of the cosmological parameters 
(see also Ref. 40 ). 

5 Results and Discussion 

We first perform string simulations with a string sampling comoving spacing of 
1/1000 of the simulation box sizes. The dynamic ranges cover from 0.05 to 300 ?7 cq , 
with each single run having a dynamic range of at least 10 in conformal time. We 
then perform the structure formation simulations with comoving box sizes rang- 
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Fig. 13. The accuracy of the rescaling scheme (55) for cosmic string models. S ap = Ank^V (k) 
is obtained from a flat A = model using (55), while SVmm is the numerical result by solving 
perturbation equations with various choices of fi m o and SIao ■ 
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Fig. 14. The relative contribution in the power spectrum from strings at different epochs: the 
total (solid line), the early times (dot-dashed line), around r] cq (dashed line), and the late times 
(dotted line). The numbers in brackets specify in rj eq the initial and final time of each epoch. 

ing from 5-130/i _1 Mpc today, and a resolution of 128 3 -512 3 . Our simulations were 
carried out on the UK Computational Cosmology Consortium supercomputer COS- 
MOS, a Silicon Graphics Origin2000 with 20 Gbytes main memory. The code was 
parallelized to enhance performance. We also used the SGI math libraries to im- 
plement the Fast Fourier Transform. 

5.1 The semi-analytic fit 

In Fig. 15 we plot the CMBR normalized linear power spectrum (G^ 6 = G\i x 10 6 = 
1.7, 12 the most recent COBE normalization for strings) induced by long strings in 
an f2 c0 = 1 CDM cosmology with h = 0.7. The central set of numerical points was 
sourced by string network simulations beginning at 77 = 0.4?7 oq which were continued 
for 1318 expansion times in simulation boxes ranging from 32-T28fo~ 1 Mpc, with a 
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Fig. 15. The comparison of our semi-analytical fit and our simulation result for the CDM model 
excluding loops. The top-right, central, and bottom-left solid lines with crosses are the simulation 
results in the deep radiation, transition through to deep into the matter, and deep matter eras 
respectively. The dashed lines are our semi-analytical fits corresponding to the same dynamical 
ranges of those simulations. They show a good agreement with each other. The solid line is the 
semi-analytical fit with a full dynamic range from r}\ = to today. 



maximum resolution of 512 3 grid points. Given the dynamic range limitations we 
have also plotted the semi-analytic fit (47) over the full range of wave-numbers. 
The good agreement with the semi-analytic model is illustrated by the dashed line 
fits to central points, as well as to the short normalization runs in the matter and 
radiation eras (also shown in Fig. 15). We can see the little discrepancy between 
the semi-analytic fits and simulations on both the large-scale and the small-scale 
ends. This resulted from the fact that towards the beginning and the end of the 
simulations, the UETC does not fully contribute to the simulation power spectrum 
(see eq. (49)). Nevertheless, at large scales (k ~ 0.03) the full-dynamic-range run 
of the semi-analytic fit (the solid line) is well constrained by the deep-matter-era 
simulation. Given this close correspondence, which was also exhibited in HDM 
simulations, we have considerable confidence that this approximation can reproduce 
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the correct shape and amplitude of the string simulation power spectrum. 

Fig. 16 shows how the radiation perturbations affect the matter perturbation 
power spectrum. The solid line is obtained using eqns. (1) and (2) with A = K = 0, 
while the dashed line is obtained by setting S T = in cqn. (1) and hence ignoring 
eqn. (2). We can see that the radiation perturbations have effect only on inter- 
mediate scales. They boost the power spectrum by about 30% at most on scales 
k <~ 0.2/iMpc _ . This is because at very small scales, radiation will oscillate many 
times per expansion time and will have little net effect on the matter, while at very 
large scales, most power is contributed from the matter era in which the effect from 
radiation is negligible. Although ignoring 6 T has its convenience for numerical pur- 
poses, the assumption of a smooth radiation background would have ignored some 
intermediate-scale power expected to result from radiation perturbations during 
radiation-matter transition. 




k (h/Mpc) 

Fig. 16. The effect of radiation perturbations in the cosmic string-seeded CDM spectrum. The 
dynamic range goes from n = to today. 
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The results presented here are obtained by integrating (47), with !F{k,ri) pro- 
vided by (54) and £ c (fc; r/o, ??) numerically obtained from (1) and (2). 

5.2 The O m0 = 1 and A = model 

Consider first the Sl c o = 1 CDM model (I). An empirical formulae which can re- 
produce our result excluding loops within a maximum of 10% error for k = 0.01- 
100/iMpc -1 is obtained as: 

ScDM(fc) - (G M6 ) 2 41.6(0.7 g )^ , (58) 

where q — k/T and 

^ = 3 - 9 - l + (2 2 8 7 g) -o. 44 - ^ 
This fit is obtained using h = 1. The T here is called the 'shape parameter', which 
accounts for the rescaling in k due to different choices of cosmological parameters. 
In the Qbo = cosmologies, T is simply fl m oh, which reflects the scaling in the 
horizon size at Tfeq- As for the fieo 7^ models, the shape parameter T will be also 
a function of Obo as we shall discuss later. 

We calculated the standard deviation as by convolving our CDM perturbations 
with a spherical window of radius 8/i _1 Mpc to find cr 8 ( s i m )(/i = 0.5) = 0.32G^ 6 , 
cr 8(sim)(^ = 0.7) = 0.39G/X6 and cr 8 ( sim )(ft. = 1.0) = 0.47G^6- A comparison with 
the observational data points is shown at the top panel of Fig. 17. It indicates that 
strings appear to induce an excess of small-scale power and a shortage of large- 
scale power, that is, the O c o = 1 string model excluding loops requires a strongly 
scale-dependent bias. This is not necessarily a fatal flaw for the model on small 
scales because, as the corresponding HDM spectrum indicates, such excess power 
can be readily eliminated in a mixed dark matter model. However, the problem 
is less tractable on large scales where biases up to crioofobsj/^iootsim) ~ 3.9 around 
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Fig. 17. The spectra of cosmic strings (excluding loops) with CDM (left) and HDM (right) models 
in different background cosmologies. The solid lines are the fiat models with n m o + Qao = 1; the 
dashed lines are the open models with Hao = 0. In an open universe with A = 0, the spectrum has 
more large-scale power and less small-scale power as Q m o decreases. In a flat universe made by A, 
the spectrum has a scale-independent boost over the previous case with the same U m o- Here we 
use h = 0.7 and G/i6 = 1.7, 12 which is the most recent COBE normalization for cosmic strings. 
The data points with error bars are the reconstructed linear spectrum by Peacock and Dodds. 41 
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100/i~ 1 Mpc might be inferred from the data points (using G/x 6 = 1.7 and h = 0.7). 
Should we, therefore, rule out string models on this basis? a While acknowledging 
that the model (I) linear spectrum looks unattractive, there are still three important 
mitigating factors. First, the determination of the power spectrum on large scales 
around 100/i _1 Mpc remains uncertain (e.g. the observational data in Fig. 17 has 
at least a 40% normalization uncertainty). Even with the normalization based on 
new galaxy surveys 42,43 and cluster abundance 44 , this uncertainty still remains as 
large as 40% at 95% confidence level 45 . Nevertheless, it will be superseded by 
much larger more reliable data sets in the near future. Secondly, the immediate 
nonlinearity of string wakes indicates that strong biasing mechanisms might operate, 
as illustrated on large scales in the post-recombination hydrodynamic simulations 
of Rcf. 46 . Thirdly, cosmic string loops have not been included here, and we will see 
how they release this problem later. Finally, unlike inflation, defect models have 
never been wedded to an Q = 1 cosmology. 

As for the HDM results, the small-scale power (k > 0.5/iMpc -1 ) is removed by 
the neutrino free-streaming, while on large scales it remains the same as that of 
the CDM model. An empirical formulae (excluding the effect of loops and baryons) 
within a maximum of 10% error for k = 0.01-lOOft.Mpc -1 is given by: 

708 

Shdm (k) = (^ 6 ) 2 [(0 7(z) _, 2 + 5 3#3 , (60) 

where again q = k/T and T is the shape parameter. This fit again has used h — 1. 
Fig. 18 shows a comparison between the mass distributions of CDM and HDM 
models. They are slices of a simulation box of size (128/i _1 Mpc) 3 , with exactly 
the same string source. They look the same on larger scales, but the HDM model 

a In Ref. 10 , the authors used a numerical fit for the observational data to obtain Cioo(obs) = 
3.5 X 10~ 2 , and compared this with their simulated spectrum to yield bioo = 5.4. They concluded 
that cosmic string theory is ruled out because of this high bias factor. However, according to the 
observational data 41 (see Fig. 17), there arc only a few data points (with big error bars) below 
the scale of 100/i -1 Mpc, so this seems a potentially unreliable way to estimate crioo(obs)- 
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has no fine structures on smaller scales due to the neutrino free-streaming. Their 
non-Gaussian features can also be seen. Ref. 7 has investigated their non-Gaussian 
properties in detail, and concluded that on scales smaller than l.5(Qh 2 ) 1 Mpc, 
perturbations seeded by cosmic strings are very non-Gaussian. These scales may 
still be in a linear or mildly non-linear regime in an open or A-universe with T = 
Slh < 0.2. 

5.3 The cases f2 m0 < 1 and A^O 

We can observe from Fig. 17, that for open or A-cosmologies with f2 m0 w 0.1 
0.3, the string + CDM power spectrum is much more encouraging. We find that 
the bias on large scales is much less scale-dependent. For example, in model (IV) 
(O m o = 0.15, f^Ao = 0.85), the relative bias remains <xioo(obs) /^loofsim) ~ 1.4±0.2 at 
100/i _1 Mpc. In Fig. 19, we plot a% for the full gamut of open and A models. o^sim) 
induced from our simulations with GfiQ normalized by COBE 12 ' 9 is presented as 
the dot-dashed (h = 1.0), the solid (h = 0.7) and the dashed (h = 0.5) lines; CT8(obs) 
is the shaded area (95% confidence region), which is recommended by Wu, 45 and 
which incorporates the latest techniques (e.g., the inclusion of non-spherical-collapsc 
models as an improved alternative to the Press-Schechter formalism) and large-scale 
data (new galaxy surveys and cluster data) to improve earlier results. 47 ' 48 ' 49 We 
can see from Fig. 19 that the consistency test between cr 8 ( obs ) and o- 8 ( sim ) indicates, 
for h = 1.0,0.7,0.5 respectively, that 0.15 > O m0 > 0.35, 0.2 > O m0 > 0.6 and 
Jl m o > 0.4 when fi m0 + fl A0 = 1, and that 0.2 > fl m0 > 0.45, 0.3 > n m0 > 0.7 
and fimo > 0.45 when Qao = 0. For h = 0.7 in both A and open models, the 
ratio o'8(obs)/ cr 8(sim) < 2 for all f2 m o > 0.1. Combined with an analysis similar to 
Fig. 17, we found that for T = Cl m oh = 0.1-0.2, both 08( s im) and the shape of 
spectrum induced from cosmic string model match observations within acceptable 
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Fig. 18. A comparison between CDM (top) and HDM (bottom) results. They are taken from 
slices in simulation box of size(128h _1 Mpc) 3 , with exactly the same string source. They look the 
same on large scales, but the HDM model has no fine structure. The color scheme is chosen so 
that the non-Gaussian feature can be clearly distinguished. 
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uncertainties. Hence, an open or A-cosmology in the context of string + CDM 
model seems to show remarkable agreement between observations and the results 
of our numerical simulations. Indeed, similar comparisons made in an inflationary 
context with models like I-V appear to require a more strongly scale-dependent 
bias. 36 
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Fig. 19. The comparison of the observationally inferred standard deviation at the scale 8 h~ 1 Mpc, 
°"8(obs)i an d that induced from our simulation, o~$f s i m ). We have used the COBE normalization 
G/Je(r2 c o = 1,Oao = 0) = 1.7. 12 ' 9 o" g ( obs ) is shown as the shaded area; 45 o" g ( sim ) is plotted as 
dot-dashed (h = 1.0), solid (h = 0.7) and dashed(/i = 0.5) lines. 

As for the HDM results, the comparison with observation seems to require a 
strongly scale-dependent bias for any choice of the cosmological parameters (models 
I-V). However, the lack of small-scale power may be partially overcome, at least, 
if baryons are properly included in the analysis. Further investigation using a 
hydrodynamical code will be required to investigate if galaxies form early enough. 

5.4 Loop effect 

Now let us include the effect from loops. Figure 20 shows the power spectrum of 
density perturbations induced by long strings and by cosmic string loops for / = 1 
with a small dynamic range from 2.5 to 5ry eq . We can see that when compared 
with the spectrum induced by static loops (dot dashed), the amplitude of small- 
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Fig. 20. Small dynamic range power spectra of density perturbations seeded by long strings (thick 
solid), by loops with initial velocities v* switched to zero (dot-dashed), and by loops with u* 
determined by string network evolution (dashed and thin solid). The thin solid line includes the 
effect of gravitational decay of the loop energy, while the other two loop lines don't but with loops 
removed after a period of time r* = t* . 

scale perturbations induced by moving loops (dashed) is clearly reduced by their 
motion. However, their large-scale power is higher because of the dependence of the 
gravitational interaction on the loop velocities, especially when they are relativistic. 
We also see that the gravitational decay of loop energy (thin solid) damps the overall 
amplitude of the power spectrum (dashed) by about a factor of 3. We notice that 
between the long-string correlation scale w 20/t] (see subsection 3.2) and the 
scale /cl ~ 10fc|, the slope of the loop spectrum (thin solid) is exactly the same as 
that of the long-string spectrum n w —2.25. We have shown in subsection 3.3 that 
this close correspondence is due to copious loop production being strongly correlated 
with long string intercommuting events and the collapse of highly curved long-string 
regions, 14 that is, near the strongest long-string perturbations. Moreover, these 
correlations persist in time with the subsequent motion of loops and long strings 
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lying preferentially in the same directions, a phenomenon which has been verified by 
observing animations of string network evolution. These correlations between loops 
and long strings, however, have a lower cutoff represented by the mean loop spacing 
g?l ~ 1 . Below cJl, the effects of individual filaments swept out by moving loops 
can be identified. In terms of the power spectrum, for k < k^ the loops are strongly 
correlated with the long strings and therefore reinforce the wake-like perturbations, 
while for k > k^ their filamentary perturbations increase the spectral index by 
about one tons —1.25; this change is expected on geometrical grounds. 

In figure 21 we plot the power spectra of density perturbations seeded by long 
strings Voc(k), by small loops Vh(k), and by both loops and long strings Vt t(k). 
The dynamic range here extends from 0.6 to 7.5rj eq . As expected Ph(k) scales more 




Fig. 21. The lower set of 3 lines are Vj J (k) for / = 0.5 (dotted), 1 (dot-dashed) and 2 (dashed). 
"Poo(fc) is plotted as a solid line. The upper set of lines are Vtot{k) with corresponding line styles 
and / values to the lower set of lines. 



moderately than f 2 but more strongly than / (see (37)). It is also apparent that 
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the perturbations induced by long strings and by loops are positively correlated 
with Ptot(k) > Vhik) + Voo(k) throughout the whole scale range. This positive 
correlation between loops and long strings boosts the large-scale Voo(k) by a factor 
of 1.5, 1.8 and 2.2 to reach Vtot(k) for / = 0.5, 1 and 2 respectively, even if Vh(k) 
is a relatively small fraction of Voo{k) on these scales. 

Figure 22 shows the correlation coefficient between the long-string and loop 
induced perturbations: 

(5 2 o) l/ 2((5 2 )1/2 • ^ 

We see that long strings and loops are strongly positively correlated on large scales, 
but weakly correlated on small scales where the loops dominate the perturbations 
(see figure 21). The threshold k t between these two regimes must be significantly 
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Fig. 22. The correlation coefficient between the long-string and loop induced perturbations, with 
/ = 0.5, 1, 2, 4, 6 (downwards). 



larger than fc L because, for k < fc L , V\,{k) is well below and roughly parallel to 
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'Poo(k) (see figures 20 and 21). We also verify that Vtot{k) /Poo(k) is approximately 
a constant for k < <k t , which again provides strong evidence for the fact that 
loops behave as part of the long-string network on large scales. 

Given these properties of the string power spectra, one can easily construct a 
semi- analytic model for Vtot(k) as for V 00 (k). 6 ' 8 We first multiply the structure 
function T(k,r]) of Voo(k) by J{r},f) to account for the boost Ptot(fc) /'Poo(fc) on 
large scales (k < k t ): 



where J(f) = 0.806, 1, 1.22, 1.44 and 1.56 for / = 0.5, 1, 2, 4 and 6 respectively. 
We then multiply it again by a numerically verified form 



to account for the turnover for k > fc t (r/, /). J(i],f) 1S calibrated phenomeno- 
logically from simulations deep in the radiation era through to those deep in the 
matter era. In the pseudo-scaling regime for the loop size, fc t is revealed to be at 
least 10fc{ ~ 2OO/77 depending on /. Thus we can carry out a full-dynamic-range 
integration to obtain Vtot{k). In figure 23 we compare this Vt t(k) and r P oc (k) with 
observations. 6 ' 8 The background cosmology is Q c = 0.15, f^A = 0.85 and h = 0.7, 
and we have used the COBE normalization G/j = 1.7 x 10~ 6 throughout. 12 Since 
loops are point-like and they have little impact through the Kaiser-Stcbbins effect 
on COBE-scale CMBR anisotropies, we expect this normalization to be very weakly 
dependent on the value of /; indeed, loops were found to produce a negligible ef- 
fect in Ref. 50 . Thus we see from figure 23 that for / > 0.5, loops can contribute 
significantly to the total power spectrum and ease the large-scale bias problem seen 
previously. 6 ' 8,10 Definite conclusions, therefore, about biasing in cosmic string mod- 
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Fig. 23. Comparison of the observational power spectrum 41 with Vaa(k) (solid), and "Ptot(fc) for 
/ = 0.5, 1, 2, 4, 6 (dashed, upwards), with a full dynamic range. 



els will need further advances in determining the magnitude of the parameter /, 
while all future large-scale structure simulations will now require the the inclusion 
of loops. 

These additional complications in modeling cosmic string structure formation 
arc most obvious on small scales, where even higher resolution and large dynamic 
range simulations will be required. Within the present pseudo-scaling regime for 
loop size, wc know that k t > k L > lOfc^ as shown in figure 20 and discussed 
previously. Taking this extreme minimum k t = then, wc find that the semi- 

analytic model over the full dynamic range gives at most a 2% difference in Vtot(k) 
for k < 1/iMpc -1 when the filament term H{k,-q,f) is excluded from T{k,r]) (for 
il c = 0.15, £1a = 0.85 and h = 0.7). This means that although the simulations 
described in this letter are already on the verge of present computer capabilities, a 
further detailed study on small scales will improve only the overall normalization 
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of Vtot(k) but not the shape revealed here if / is constant throughout. We notice 
that although / should be constant in both deep radiation and deep matter eras 
due to the scaling of string network, it may vary over the transition era. This time 
dependence in / could slightly alter the shape of Vtot(k) given here, and this deserves 
further investigation on cosmic string evolution. We also note that advances in 
understanding loop formation mechanisms will also be crucial in quantifying the 
importance of the gravitational radiation background emitted by a cosmic string 
network and its effect on large-scale structure and CMBR anisotropics. 51 

5.5 CMBR anisotropies 

A key feature of all these string-induced power spectra is the influence of the slow 
relaxation to the matter era string density from the much higher radiation string 
density, which has an effective structure function !F(k,r]) in (54) with approxi- 
mately 2.5 times more power than the matter era version. Even by recombination 
in an f2 m0 = 1 (h = 0.7) cosmology, the string density ((i]) is more than twice its 
asymptotic matter era value to which we normalize on COBE scales (see Fig. 8). 
This implies that the string model provides higher than expected large-scale power 
around 100/i _1 Mpc and below. Interestingly, this can also be expected to produce a 
significant Doppler-like peak on small angle CMBR scales, an effect noted in Ref. 12 
but not observed because matter era strings were employed. In Ref. 5 , global strings 
were evolved through the transition era without a strong Doppler peak emerging; 
however, this depended on a field theory simulation of limited dynamic range. Since 
global strings on cosmological scales behave more like local strings, the linear power 
spectrum approach presented here potentially should more accurately represent the 
actual global string power spectrum. Recent work in Ref. 23 using a phcnomcno- 
logical semi-analytic approach confirms that such Doppler-like features can result 
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from significant non-scaling effects during the transition era. Especially in the light 
of seeing the first Doppler peak in the recent observations from MAXIMA- 1 52 and 
BOOMERANG 53 , the cosmic-string-predicted first peak requires more careful and 
detailed study. We should also note that the currently observed CMB power spec- 
trum has quite a narrow first peak and no evident secondary peaks 54 , and this may 
imply a hybrid scenario of structure formation, which is a mix of both inflation and 
topological defects 55 . 

5.6 Uncertainties 

Finally, it is appropriate to comment on the key uncertainties affecting these cal- 
culations. To summarize at the outset, these uncertainties primarily affect the am- 
plitude of the string power spectrum, rather than its overall shape which is a more 
robust feature. The most recent COBE string normalization is G^lq w 1.7, 12 which 
is at variance with a previous COBE normalization GfiQ s=s 1 using the same string 
simulation. 50 This remains to be satisfactorily resolved at the time of writing, but 
systematic relativistic effects seem likely to have affected the earlier result. Next 
there is uncertainty implementing the effective compensation scale at which the 
perturbations are cut-off on large scales, which we have discussed already. Thirdly, 
there are uncertainties in the long-string energy density as we have addressed in sec- 
tion 3.2. This causes an uncertainty of about 10% in the resulting power spectrum. 
All these uncertainties add up to a factor of about 2 in the final power spectrum 
while remaining the overall shape largely unchanged. Finally, the high correlation 
between loops and long strings always boosts the power spectrum up to a factor 
of 3, if the loop lifetime is not much smaller than the Hubble time. Therefore, to 
include the effect from loops has become necessary for any further development in 
studying structure formation by cosmic string models. 
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6 Conclusion 

In this paper we have described the results of high-resolution numerical simula- 
tions of structure formation seeded by a local cosmic string network with a large 
dynamical range reporting at length, for the first time, on the effect of loops and 
modifications due to the radiation-matter transition. In the regime of large-scale 
structure formation, the most serious problem with both cosmic strings and the 
standard CDM model (inflation) is that they produce too much small-scale power 
and insufficient large-scale power. From the results and discussion in section 5.3, 
we see for cosmic strings this problem can be relaxed a great deal by considering 
an open universe or the existence of a small cosmological constant. There is an- 
other alternative which can improve the situation. Mixed dark matter scenarios 
can be employed by using (44) to smear out part of the small-scale perturbations, 
although we see the pure HDM power spectrum requires a strongly scale-dependent 
bias either on small or large scales. 

As to the shortage of the overall amplitude in the string power spectrum when 
compared with observations, we have seen that this can be overcome by including 
cosmic string loops with a lifetime comparable to the Hubble time or greater. We 
have shown that on large scales the loops behave like part of the long-string network 
and can therefore contribute significantly to the total power spectrum of density 
perturbations. At present, the typical size and lifetime of loops formed by a string 
network remains to be studied in more detail; the problem is both computationally 
and analytically challenging. However, within the scale range of interest further 
developments in this area have the potential to affect the overall amplitude of the 
spectrum, while leaving the shape largely unchanged. 

We conclude that although more work needs to be done, notably in improving 
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the implementation of compensation and in the study of loops, the picture which 
emerges for the large-scale structure power spectrum is encouraging. 
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Appendix A. Conventions and Cosmological Background Dynamics 

We assume that the universe is homogeneous and isotropic, and is filled with two 
fluids, radiation and dark matter, whose stress-energy tensors are also homoge- 
neous and isotropic on average. We will use subscripts m, c, h, r to denote dart 
matter, cold dark matter (CDM), hot dark matter (HDM), and radiation respec- 
tively, a subscript "eq" to denote the epoch of radiation-matter density equality, 
and a subscript "0" to denote the epoch today. We will ignore the contribution of 
the defect field stress energy, which is always much smaller than the total energy 
density of radiation and matter. We use Greek letters as space-time indices (e.g. 
p = 0,1,2,3), and Latin letters as spatial indices (e.g. i — 1,2,3). The metric 
signature is ( — h ++) and units used are normalized to K = c = ks = 1- 

In a flat Friedmann- Robertson- Walker (FRW) universe with only radiation and 
matter components which evolve independently and adiabatically, the scale factor 
a(r]) is determined by the unperturbed Einstein equation 

• 2 , t>- 2 87rG/? m oao A 4 

ar + ifa = (l + a) + — a 1 , (A.l) 

where a dot represents a derivative with respect to the conformal time 77, K is the 
curvature, p m is the matter energy density, A is the cosmological constant, and we 
have normalized a cq = 1. If we define f2 m = 8irGp m /3H 2 , fl r = 8nGp r /3H 2 = 
8TrGp m /3aH 2 , Q A = A/3i? 2 , and Q K = -K/a 2 H 2 , where H = a/a 2 is the Hubble 
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parameter, then we have from (A.l) that fi m + fl r + Cl\ + = 1 and 

^ao _ A_ = -3K 

^mo 87rGp m0 ' fi m o 87rGp m0 aQ 

We notice that f2 r0 /f2 m o = a 1 < 1. We also define 

^KO ^AO 



(A.2) 



C 



(A.3) 



where B, C <C 1 according to the current observational results. Thus we can rewrite 
(A.l) as 



da \ ' 
d^J 



A 2 (l + a + Ba 2 +Ca 4 ), 



where 



i./ 1 

^eq Jo (1 + « 



da 



A. 



(A.4 



(A.5) 



VcqJo U + a + Bo? + Ca^) 1 / 2 
(A.4) can then be numerically evaluated with certain choice of f2 m o, ^ao and 
Hko- Assuming three species of neutrinos and because at ?7 oq both the curva- 
ture and the cosmological constant terms are negligible in (A.l), we obtain clq = 
23219fi m0 /i 2 , r/ oq = 16.3098 (f2 m o/i 2 ) _1 Mpc today, and the physical time t cq — 
3.4058 x 10 10 (f2 m o/i 2 ) -2 scc. In certain cases, (A.4) can be exactly solved: 

(i) K = A = (i.e. a„o = 1, O A0 = 0): 



a(r]) = A 2 r] 2 /4 + Ar) , 
t{rj) = AV/12 + V/2, 



(A.6) 
(A.7) 



which give ?7 cq = 3i cq / V2- 
(ii) K < 0, A = (i.e. fi m0 < 1, £l A0 = 0): 

a(f]) = [cosh {AVBri) + 2\f~B sinh (AVBrj) - 1 



t(v) = 



2B 
1 

AB 



cosh (Ay/Brj) + ^= sinh (A\/Brj) - ^ - 1 



(A.8) 
(A.9) 
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Figure A.l: The effect of the curvature and the cosmological constant in the back- 
ground cosmology. Plotted are the exact solutions of scale factor a{-q) (see text 
for details). The square, triangle, circle and diamond mark the epoch today (with 
h = 0.7) for different models. 



(hi) K > 0, A = (i.e. Q m0 > 1, ^ A0 = 0): 



™ = Jb 



cos (AV^Btj) - 2v / ^Bsin (AV^Brj) - 1 
1 



cos (Av — Br/) 



At] 



. u >»> {Av—Br}) ^- — 1 



(A.10) 
(A.ll) 



We notice that at early times (A. 8, A. 9) and (A. 10, A. 11) decay to (A. 6, A. 7) exactly. 
At late times (A. 6, A. 7), (A. 8, A. 9) and (A. 10, A. 11) give the asymptotic forms a cx 
rj 2 , exp(A\/Bri) and 1 — cos( J 4\/— Brf) (before the recollapse), or a a a cx t 
and a cx 1 — cos[2A(— B) 3 / 2 t] (before the recollapse) respectively. Figure A.l shows 
some examples of these solutions. 
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Appendix B. Power Spectrum Conventions 



Wc define the Fourier transform as: 



/(k) = I J d 3 x/(x)e 



(B.l) 



where k is the wave vector, x the physical coordinates, and the integration taken 
over a large volume V. We will use a tilde ~ to denote the Fourier transform 
of a function. The power spectrum V(k) for a mass density distribution p(x) is 
then defined as the spherically symmetric Fourier transform of the autocorrelation 
function £(r) of density contrast 5(x) = 5p/p av : 



where k = |k|, r is the correlation distance, and 5^ is the Fourier transform of 
the perturbations <5(x). A useful quantity when comparing numerical results with 
observation is the dimensionless "spectrum of the matter density contrast" , which 
is given by contributions of Sk over a logarithmic interval Ak/k <~ 1: 



This quantity is dimensionless and gives the rms density fluctuation on a particular 
length-scale I = 2n/k. Another useful statistic is the mass fluctuation amplitude at 
certain length-scale A, i.e. the standard deviation of the mass density distribution 
<j\. We convolve the density contrast <5(x) with a spherical window of radius A to 
obtain the variance: 



a\ = ^\ \w(k\)\ 2 k 2 V(k)dk with w(x) = Hs[nX 3 XC ° SX) . (B.4) 



We note that w(x) is the Fourier transform of a spherical window and V(k) in (B.4) 
can be calculated from (B.2). 




(B.2) 




(B.3) 



